Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing.

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Persistent SARS-CoV-2 infection with accumulation of mutations in a patient with poorly controlled HIV infection.

A 22-year-old female with uncontrolled advanced HIV infection was persistently infected with SARS-CoV-2 beta variant for 9 months, the virus accumulating >20 additional mutations. Antiretroviral therapy suppressed HIV and cleared SARS-CoV-2 within 6-9 weeks. Increased vigilance is warranted to benefit affected individuals and prevent the emergence of novel SARS-CoV-2 variants.

SARS-CoV-2 Genetic Diversity and Lineage Dynamics in Egypt during the First 18 Months of the Pandemic.

COVID-19 was first diagnosed in Egypt on 14 February 2020. By the end of November 2021, over 333,840 cases and 18,832 deaths had been reported. As part of the national genomic surveillance, 1027 SARS-CoV-2 near whole-genomes were generated and published by the end of July 2021. Here we describe the genomic epidemiology of SARS-CoV-2 in Egypt over this period using a subset of 976 high-quality Egyptian genomes analyzed together with a representative set of global sequences within a phylogenetic framework. A single lineage, C.36, introduced early in the pandemic was responsible for most of the cases in Egypt. Furthermore, to remain dominant in the face of mounting immunity from previous infections and vaccinations, this lineage acquired several mutations known to confer an adaptive advantage. These results highlight the value of continuous genomic surveillance in regions where VOCs are not predominant and the need for enforcement of public health measures to prevent expansion of the existing lineages.

Omicron BA.4/BA.5 escape neutralizing immunity elicited by BA.1 infection.

SARS-CoV-2 Omicron (B.1.1.529) BA.4 and BA.5 sub-lineages, first detected in South Africa, have changes relative to Omicron BA.1 including substitutions in the spike receptor binding domain. Here we isolated live BA.4 and BA.5 viruses and measured BA.4/BA.5 neutralization elicited by BA.1 infection either in the absence or presence of previous vaccination as well as from vaccination without BA.1 infection. In BA.1-infected unvaccinated individuals, neutralization relative to BA.1 declines 7.6-fold for BA.4 and 7.5-fold for BA.5. In vaccinated individuals with subsequent BA.1 infection, neutralization relative to BA.1 decreases 3.2-fold for BA.4 and 2.6-fold for BA.5. The fold-drop versus ancestral virus neutralization in this group is 4.0-fold for BA.1, 12.9-fold for BA.4, and 10.3-fold for BA.5. In contrast, BA.4/BA.5 escape is similar to BA.1 in the absence of BA.1 elicited immunity: fold-drop relative to ancestral virus neutralization is 19.8-fold for BA.1, 19.6-fold for BA.4, and 20.9-fold for BA.5. These results show considerable escape of BA.4/BA.5 from BA.1 elicited immunity which is moderated with vaccination and may indicate that BA.4/BA.5 may have the strongest selective advantage in evading neutralization relative to BA.1 in unvaccinated, BA.1 infected individuals.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Omicron infection enhances Delta antibody immunity in vaccinated persons.

The extent to which Omicron infection1-9, with or without previous vaccination, elicits protection against the previously dominant Delta (B.1.617.2) variant is unclear. Here we measured the neutralization capacity against variants of severe acute respiratory syndrome coronavirus 2 in 39 individuals in South Africa infected with the Omicron sublineage BA.1 starting at a median of 6 (interquartile range 3-9) days post symptom onset and continuing until last follow-up sample available, a median of 23 (interquartile range 19-27) days post symptoms to allow BA.1-elicited neutralizing immunity time to develop. Fifteen participants were vaccinated with Pfizer's BNT162b2 or Johnson & Johnson's Ad26.CoV2.S and had BA.1 breakthrough infections, and 24 were unvaccinated. BA.1 neutralization increased from a geometric mean 50% focus reduction neutralization test titre of 42 at enrolment to 575 at the last follow-up time point (13.6-fold) in vaccinated participants and from 46 to 272 (6.0-fold) in unvaccinated participants. Delta virus neutralization also increased, from 192 to 1,091 (5.7-fold) in vaccinated participants and from 28 to 91 (3.0-fold) in unvaccinated participants. At the last time point, unvaccinated individuals infected with BA.1 had low absolute levels of neutralization for the non-BA.1 viruses and 2.2-fold lower BA.1 neutralization, 12.0-fold lower Delta neutralization, 9.6-fold lower Beta variant neutralization, 17.9-fold lower ancestral virus neutralization and 4.8-fold lower Omicron sublineage BA.2 neutralization relative to vaccinated individuals infected with BA.1. These results indicate that hybrid immunity formed by vaccination and Omicron BA.1 infection should be protective against Delta and other variants. By contrast, infection with Omicron BA.1 alone offers limited cross-protection despite moderate enhancement.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Targeted Sanger sequencing to recover key mutations in SARS-CoV-2 variant genome assemblies produced by next-generation sequencing.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is adaptively evolving to ensure its persistence within human hosts. It is therefore necessary to continuously monitor the emergence and prevalence of novel variants that arise. Importantly, some mutations have been associated with both molecular diagnostic failures and reduced or abrogated next-generation sequencing (NGS) read coverage in some genomic regions. Such impacts are particularly problematic when they occur in genomic regions such as those that encode the spike (S) protein, which are crucial for identifying and tracking the prevalence and dissemination dynamics of concerning viral variants. Targeted Sanger sequencing presents a fast and cost-effective means to accurately extend the coverage of whole-genome sequences. We designed a custom set of primers to amplify a 401 bp segment of the receptor-binding domain (RBD) (between positions 22698 and 23098 relative to the Wuhan-Hu-1 reference). We then designed a Sanger sequencing wet-laboratory protocol. We applied the primer set and wet-laboratory protocol to sequence 222 samples that were missing positions with key mutations K417N, E484K, and N501Y due to poor coverage after NGS sequencing. Finally, we developed SeqPatcher, a Python-based computational tool to analyse the trace files yielded by Sanger sequencing to generate consensus sequences, or take preanalysed consensus sequences in fasta format, and merge them with their corresponding whole-genome assemblies. We successfully sequenced 153 samples of 222 (69 %) using Sanger sequencing and confirmed the occurrence of key beta variant mutations (K417N, E484K, N501Y) in the S genes of 142 of 153 (93 %) samples. Additionally, one sample had the Y508F mutation and four samples the S477N. Samples with RT-PCR Ct scores ranging from 13.85 to 37.47 (mean=25.70) could be Sanger sequenced efficiently. These results show that our method and pipeline can be used to improve the quality of whole-genome assemblies produced using NGS and can be used with any pairs of the most used NGS and Sanger sequencing platforms.

SARS-CoV-2 prolonged infection during advanced HIV disease evolves extensive immune escape.

Characterizing SARS-CoV-2 evolution in specific geographies may help predict properties of the variants that come from these regions. We mapped neutralization of a SARS-CoV-2 strain that evolved over 6 months from ancestral virus in a person with advanced HIV disease in South Africa; this person was infected prior to emergence of the Beta and Delta variants. We longitudinally tracked the evolved virus and tested it against self-plasma and convalescent plasma from ancestral, Beta, and Delta infections. Early virus was similar to ancestral, but it evolved a multitude of mutations found in Omicron and other variants. It showed substantial but incomplete Pfizer BNT162b2 escape, weak neutralization by self-plasma, and despite pre-dating Delta, it also showed extensive escape of Delta infection-elicited neutralization. This example is consistent with the notion that SARS-CoV-2 evolving in individual immune-compromised hosts, including those with advanced HIV disease, may gain immune escape of vaccines and enhanced escape of Delta immunity, and this has implications for vaccine breakthrough and reinfections.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa.

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Omicron extensively but incompletely escapes Pfizer BNT162b2 neutralization.

The emergence of the SARS-CoV-2 variant of concern Omicron (Pango lineage B.1.1.529), first identified in Botswana and South Africa, may compromise vaccine effectiveness and lead to re-infections1. Here we investigated Omicron escape from neutralization by antibodies from South African individuals vaccinated with Pfizer BNT162b2. We used blood samples taken soon after vaccination from individuals who were vaccinated and previously infected with SARS-CoV-2 or vaccinated with no evidence of previous infection. We isolated and sequence-confirmed live Omicron virus from an infected person and observed that Omicron requires the angiotensin-converting enzyme 2 (ACE2) receptor to infect cells. We compared plasma neutralization of Omicron relative to an ancestral SARS-CoV-2 strain and found that neutralization of ancestral virus was much higher in infected and vaccinated individuals compared with the vaccinated-only participants. However, both groups showed a 22-fold reduction in vaccine-elicited neutralization by the Omicron variant. Participants who were vaccinated and had previously been infected exhibited residual neutralization of Omicron similar to the level of neutralization of the ancestral virus observed in the vaccination-only group. These data support the notion that reasonable protection against Omicron may be maintained using vaccination approaches.

Cost-effectiveness of public health strategies for COVID-19 epidemic control in South Africa

Background Healthcare resource constraints in low and middle-income countries necessitate selection of cost-effective public health interventions to address COVID-19. Methods We developed a dynamic COVID-19 microsimulation model to evaluate clinical and economic outcomes and cost-effectiveness of epidemic control strategies in KwaZulu-Natal, South Africa. Interventions assessed were Healthcare Testing (HT), where diagnostic testing is performed only for those presenting to healthcare centres;Contact Tracing (CT) in households of cases;Isolation Centres (IC), for cases not requiring hospitalisation;community health worker-led Mass Symptom Screening and diagnostic testing for symptomatic individuals (MS);and Quarantine Centres (QC), for contacts who test negative. Given uncertainties about epidemic dynamics in South Africa, we evaluated two main epidemic scenarios over 360 days, with effective reproduction numbers (R e ) of 1.5 and 1.2. We compared HT, HT+CT, HT+CT+IC, HT+CT+IC+MS, HT+CT+IC+QC, and HT+CT+IC+MS+QC, considering strategies with incremental cost-effectiveness ratio (ICER) <US$1,290/year-of-life saved (YLS) to be cost-effective. Findings With R e 1.5, HT resulted in the most COVID-19 deaths and lowest costs over 360 days. Compared with HT, HT+CT+IC+MS reduced mortality by 76%, increased costs by 16%, and was cost-effective (ICER $350/YLS). HT+CT+IC+MS+QC provided the greatest reduction in mortality, but increased costs by 95% compared with HT+CT+IC+MS and was not cost-effective (ICER $8,000/YLS). With R e 1.2, HT+CT+IC+MS was the least costly strategy, and HT+CT+IC+MS+QC was not cost-effective (ICER $294,320/YLS). Interpretation In South Africa, a strategy of household contact tracing, isolation, and mass symptom screening would substantially reduce COVID-19 mortality and be cost-effective. Adding quarantine centres for COVID-19 contacts is not cost-effective.

Clinical outcomes and cost-effectiveness of COVID-19 vaccination in South Africa.

Low- and middle-income countries are implementing COVID-19 vaccination strategies in light of varying vaccine efficacies and costs, supply shortages, and resource constraints. Here, we use a microsimulation model to evaluate clinical outcomes and cost-effectiveness of a COVID-19 vaccination program in South Africa. We varied vaccination coverage, pace, acceptance, effectiveness, and cost as well as epidemic dynamics. Providing vaccines to at least 40% of the population and prioritizing vaccine rollout prevented >9 million infections and >73,000 deaths and reduced costs due to fewer hospitalizations. Model results were most sensitive to assumptions about epidemic growth and prevalence of prior immunity to SARS-CoV-2, though the vaccination program still provided high value and decreased both deaths and health care costs across a wide range of assumptions. Vaccination program implementation factors, including prompt procurement, distribution, and rollout, are likely more influential than characteristics of the vaccine itself in maximizing public health benefits and economic efficiency.

HIV status alters disease severity and immune cell responses in Beta variant SARS-CoV-2 infection wave.

There are conflicting reports on the effects of HIV on COVID-19. Here, we analyzed disease severity and immune cell changes during and after SARS-CoV-2 infection in 236 participants from South Africa, of which 39% were people living with HIV (PLWH), during the first and second (Beta dominated) infection waves. The second wave had more PLWH requiring supplemental oxygen relative to HIV-negative participants. Higher disease severity was associated with low CD4 T cell counts and higher neutrophil to lymphocyte ratios (NLR). Yet, CD4 counts recovered and NLR stabilized after SARS-CoV-2 clearance in wave 2 infected PLWH, arguing for an interaction between SARS-CoV-2 and HIV infection leading to low CD4 and high NLR. The first infection wave, where severity in HIV negative and PLWH was similar, still showed some HIV modulation of SARS-CoV-2 immune responses. Therefore, HIV infection can synergize with the SARS-CoV-2 variant to change COVID-19 outcomes.

The emergence and ongoing convergent evolution of the SARS-CoV-2 N501Y lineages.

The independent emergence late in 2020 of the B.1.1.7, B.1.351, and P.1 lineages of SARS-CoV-2 prompted renewed concerns about the evolutionary capacity of this virus to overcome public health interventions and rising population immunity. Here, by examining patterns of synonymous and non-synonymous mutations that have accumulated in SARS-CoV-2 genomes since the pandemic began, we find that the emergence of these three "501Y lineages" coincided with a major global shift in the selective forces acting on various SARS-CoV-2 genes. Following their emergence, the adaptive evolution of 501Y lineage viruses has involved repeated selectively favored convergent mutations at 35 genome sites, mutations we refer to as the 501Y meta-signature. The ongoing convergence of viruses in many other lineages on this meta-signature suggests that it includes multiple mutation combinations capable of promoting the persistence of diverse SARS-CoV-2 lineages in the face of mounting host immune recognition.

Sixteen novel lineages of SARS-CoV-2 in South Africa.

The first severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in South Africa was identified on 5 March 2020, and by 26 March the country was in full lockdown (Oxford stringency index of 90)1. Despite the early response, by November 2020, over 785,000 people in South Africa were infected, which accounted for approximately 50% of all known African infections2. In this study, we analyzed 1,365 near whole genomes and report the identification of 16 new lineages of SARS-CoV-2 isolated between 6 March and 26 August 2020. Most of these lineages have unique mutations that have not been identified elsewhere. We also show that three lineages (B.1.1.54, B.1.1.56 and C.1) spread widely in South Africa during the first wave, comprising ~42% of all infections in the country at the time. The newly identified C lineage of SARS-CoV-2, C.1, which has 16 nucleotide mutations as compared with the original Wuhan sequence, including one amino acid change on the spike protein, D614G (ref. 3), was the most geographically widespread lineage in South Africa by the end of August 2020. An early South African-specific lineage, B.1.106, which was identified in April 2020 (ref. 4), became extinct after nosocomial outbreaks were controlled in KwaZulu-Natal Province. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify new lineages and inform measures to control the spread of SARS-CoV-2. Such genomic surveillance presented in this study has been shown to be crucial in the identification of the 501Y.V2 variant in South Africa in December 2020 (ref. 5).